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Abstract. The supersymmetric (SUSY) structure of coupled-mode equations that 
describe scattering of optical waves in one-dimensional Bragg gratings is highlighted. 
This property can find applications to the synthesis of special Bragg filters and 
distributed-feedback (DFB) optical cavities. In particular, multiple SUSY (Darboux- 
Crum) transformations can be used to synthesize DFB hlters with any desired number 
of resonances at target frequencies. As an example, we describe the design of a DFB 
structure with a set of equally-spaced resonances, i.e. a frequency comb transmission 
filter. 
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In field theory, supersymmetry (SUSY) is a model that relates fermions to bosons 
in order to obtain a unihed description of all basic interactions of nature. While 
SUSY hnds severe physical and mathematical obstacles as a realistic held theoretical 
model, SUSY quantum mechanics, i.e. the application of the SUSY superalgebra to 
quantum mechanics as opposed to quantum held theory, has ohered deeper insights 
into various aspects of standard non-relativistic quantum mechanics, creating new areas 
of research in quantum mechanics itself |T]. For example, SUSY clarihes why certain 
one-dimensional potentials are analytically solvable and others are not, helps in hnding 
new solvable potentials, exemplihes the factorisation (Darboux) method, and generates 
potentials with similar scattering properties. 

Owing to the formal analogy between the Helmholtz equation of monochromatic 
light waves in dielectric media and the stationary non-relativistic Schrodinger equation, 
the methods of SUSY have been recently transferred into the optical context [21 El IH IS] to 
synthesize dielectric optical media with tailored scattering and localization properties of 
potential interest to a wide variety of applications. SUSY optical structures can be used 
to realize global phase matching, efficient mode conversion, spatial multiplexing m , as 
well as to design transparent defects, interfaces and optical intersections HlEli. SUSY 
can be applied to discretized light in waveguide lattices as well mm- Periodic optical 
media, such as distributed feedback (DFB) structures in hbers or waveguides, provide an 
important class of optical devices. While SUSY can be applied in principle to periodic 
potentials (see, for instance m El E]), analytical results in this case are extremely 
difficult because of the lack of exact analytical forms for the energy dispersion curves and 
Bloch functions. Moreover, application of SUSY to locally-periodic potentials, in which 
the local period and potential depth undergo slow changes, is extremely challenging. 
For locally-periodic shallow gratings, light dynamics can be described with an excellent 
accuracy by couple-mode theory (GMT), where the slow changes of the grating period 
and/or grating depth are simply included in a complex scattering potential [T51 |T3]. In 
GMT, coupling of the two counter-propagating waves in the medium is described by a 
set of coupled-mode equations (GME), which can be cast in a form analogous to the 
Dirac equation of a fermionic particle subjected to an external scalar and/or vectorial 
potential (see, for instance, UniEllIH]). 

In this work we highlight the SUSY structure of GME in DFB structures, and 
show that it can be suitably exploited to synthesize Bragg grating hlters with a given 
number of resonances at desired frequencies. This task is accomplished by application 
of multiple SUSY (Darboux-Grum) transformations [T^ 120] . which lead to transparent 
potentials of Kay-Moses type [21] as a special case. Interestingly, the grating prohle 
is simply determined by the superpotential of the underlying model. As an example, 
we apply our method to the synthesis of a DFB hlter that supports discrete combs 
of transmission resonances. Application of SUSY to Bragg gratings thus provides an 
interesting analytical tool for the synthesis of DFB hlters, which is complementary 
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to other analytical synthesis methods based e.g. on the Gelf’and-Levitan-Marchenko 
inverse scattering [22l |23] . 


2. Coupled-mode equations and supersymmetry 


2.1. Coupled-mode theory 


Light propagation in nonnniform grating strnctnres, snch as those realized in optical 
fibers or wavegnides, is a well-known problem in the optics of periodic media. For 
snfficiently shallow gratings that resnlt only in small reflections within a distance of one 
wavelength, CMT is an excellent approximation to the exact problem In an 

effective one-dimensional model and for first-order Bragg scattering, the longitndinal 
refractive index modnlation can be written as n{z) = no An h{z) cospvrz/A -|- 0 ( 2 :)], 
where no is the effective mode index in absence of the grating. An <C no is a reference 
valne of the refractive index change of the grating, A is the nominal grating period 
defining the Bragg freqnency ub = 7rc/(Ano), c is the speed of light in vacnnm, 
and h{z), (j){z) describe the amplitnde and phase profiles, respectively, of the grating. 
Note that, for a non-sinnsoidal modnlation of the grating, h{x) exp[i(l){z)] is defined 
by the first-order Fonrier component of the locally-periodic refractive index [T3]. Note 
also that the CMT reprodnces with excellent accnracy the resnlts obtained by exact 
transfer matrix methods for shallow np to moderate gratings {An/no <~ 0. 1) ig, 
which is the case of corrngated DFB wavegnide filters or fiber gratings. To stndy 
Bragg scattering of connterpropagating waves at freqnencies close to ujb, let E{z,t) = 
{u{z, t) exp(— icjst -|- 2'Kizno/Xo) + v{z, t) exp(—— 27iizno/ Xq) -|- c.c.} be the electric 
field in the gniding strnctnre, where Aq = 2noA is the Bragg wavelength. The envelopes 
u and V of connterpropagating waves then satisfy the following CME mm 


i d \ d\ 

(d 1 . 
\dz-v,dtY = -^^ 


( 1 ) 

( 2 ) 


where 


q{z) = ^^h(^)exp[i0(2;)] 

is the complex scattering potential and Vg 
freqnency. 


(3) 


~ c/no is the gronp velocity at the Bragg 


2.2. Supersymmetric structure of coupled-mode equations 

It was noticed in previons works [laiiTi [T8] that Eqs.(l) and (2) can be viewed as a 
Dirac-type eqnation with scalar and vectorial conplings for the spinor {u,v), nnder a 
snitable choice of the Dirac matrices. SUSY and intertwining (Darbonx) methods for the 
Dirac eqnation have been highlighted in several works (see, for instance, m [251 [261127] 
and references therein), however their potential nse in optics of periodic media has been 
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so far overlooked (with the exception of Ref. [28].) To highlight the SUSY structure of 
CME, it is worth introducing the dimensionless variables x = zjZ and r = t/T , with 
characteristic spatial Z and time T scales dehned by 


^0 rr ^ ^ono 

- 1 = - — -. 

nAn ’ Vg wAnc 


(4) 


and the new envelopes 


'ipi{z,t) 


u{z, t) — v{z, t) 

71 


'lp2{z,t) 


u{z, t) + v{z, t) 

71 


(5) 


For monochromatic waves {u,v) ~ exp(—iUr) with a normalized frequency detuning 
E = {oj — ujb)T = 27r(z/ — i'b)T from the Bragg reference frequency ojb = 27rz/5, Eqs.(l) 
and (2) take the form 


{E — m) '^1 = iP-ip2 ( 6 ) 

{E + m)ip2 =-iP+Pi- (7) 


In Eqs.(6) and (7), the operators P± are dehned by 

+ W{x) (8) 

where m(x) and W(x) are given by 


m(x) = h{x) cos[(j){x)] , W{x) = —h(x) sin[0(x)]. (9) 

Here we focus our attention to the special class of grating prohles such that m{x) 
is a constant (i.e. independent of x), for which the SUSY structure of the Dirac 
equation follows in a rather straightforward way. The m(x)=constant case includes 
the important class of periodic gratings with no chirp, corresponding to 0(x) = —7r/2 
and thus m(x) = 0 and lU(x) = h{x) [see Eq.(9)]. The m(x) = const ^ 0 may describe, 
for example, the case of chirped gratings that realize the relativistic (Dirac) harmonic 
oscillator im. Under such an assumption, from Eqs.(6) and (7) it readily follows that 
-01 and 'ip 2 satisfy two stationary Schrodinger equations 

Hi'ijji = ei)i , ^2^2 = eV’2 (10) 


with the same eigenvalue 
e = E^ — w? 

and with partner Hamiltonians Hi and H 2 dehned by 

(f 


Hi = F_P+ = 

H 2 = P+P. = 


dx^ 

dx‘^ 


+ Vi{x) 

+ U2(x) 


( 11 ) 

( 12 ) 

(13) 


where the partner potentials Vi(x), V 2 {x) are related to the superpotential lU(x) via 
the usual relations 


/X o/ X dW , , o / X dW 

Ui(x) = lU^(x)-— , V2{x) = H^^(x) + 


dx 


dx 


( 14 ) 
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Let us now suppose that the potential Vi{x) is constructed such that the Hamiltonian 
Hi sustains a hnite number N of bound states at energies Ci, 62, 63, e^, with ei = 0 

and Cat > cn-i > ■■■ > 62 > ei. Note that, owing the factorization of Hi [Eq.(12)], the 
fundamental (lowest-order) bound state <pi{x) of Hi, Hi(j)i = ei0i = 0, is annihilated 
by the operator P+, i.e. P+<pi = 0, and the superpotential W{x) is related to 0i by the 
simple relation 


W{x) = 


{d(j)i/dx) 

01 


(15) 


The Hamiltonian H 2 is isospectral to Hi, except for the fact that at e = ei = 0 iL2 
does not sustain a bound state. Hence, apart from the e = ei = 0 case, the grating 
eigenvalue equations (6) and (7) sustain bound states at normalized frequency detunings 
= ±\/e„ + (n = 2,3, ...,N), which follows from Eq.(ll). The case e = ei = 0, 
corresponding to P = ±m, should be directly investigated by considering the original 
Eqs.(6) and (7). Since P+0i = 0, Eqs.(6) and (7) can be satisfied by taking ^jJl = 0i, 
02 = 0 and E = m. Therefore the grating sustains {2N — 1) bound states at the 
frequency detunings 


Ei = p. + rn? (n = 1, 2, 3,4,iV) 
-B” = - (n = 2,3,..,N) 


(16) 

(17) 


In particular, for a periodic grating without chirp, corresponding to 0 = —7r/2 and 
m = 0, the bound states occur at the normalized frequency detunings 


En = 0, ±-y/e2, ±\/e3, •••, Ey/e^. 


(18) 


3. Darboux-Crum transformations and synthesis of frequency comb DFB 
filters 


In this section we apply the SUSY structure of grating equations to the synthesis 
of a DFB structure that realizes a frequency comb transmission filter. For the sake 
of definiteness, we will consider the case of a periodic grating, i.e. m(x) = 0 and 
W{x) = h{x). A powerful and very general method to generate a Hamiltonian Hi 
with a desired number of bound states at target energies is the application of a 
cascade of SUSY (Darboux) transformations to a reference Hamiltonian, a procedure 
known as the Darboux-Crum theorem [T^ [20] that we briefly review here for the 
sake of completeness. Let Hq = —dj -|- Uo{x) be a reference Hamiltonian and let 
Hii%Ij{x) = P0(x), Hoipj{x) = ejipj{x) (j = 1,2, ...,Y), where ej are arbitrary numbers 
that should not necessarily belong to the point spectrum of Hq and (pj are linearly- 
independent associated functions. The functions {(pj{x),ej} are called seed functions. 
Let us then introduce the following functions 


^ yv[ipi,ip2, ...,ipN,fj]{x) 

yv[(pi,(p2, ...,(pn]{x) 

yv[(pi,(p2, ...,(pn] 


U = l,2,...,N) 


( 20 ) 








spectrum potential ^(x) h{x) spectrum potential i7j(x) h{x) 
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- 1.2 0 1-2 

normalized frequency detuning E 


Figure 1. (Color online) DFB structure that sustains 3 equally-spaced localized modes 
as obtained using the Darboux-Crum transformation with N = 2, ki = 1.2 and K 2 = 1- 
(a) Behavior of the grating amplitude h{x). (b) Behavior of the potential Ui{x). (c) 
Schematic of the grating spectrum. The three localized modes are equally spaced by 
^Ecomb = — ^2 — 0-66 and are embedded in the stop band of the grating. The 

shaded areas (|i?| > ki) correspond to allowed bands of the grating. 




(c) 


AEcom 

b 




- 

4 

-2 0 2 

4 


normalized frequency detuning E 

Figure 2. (Color online) Same as Fig.l, but for a Darboux-Crum transformation with 
iV = 3 (ki = 3, K 2 = 2.8284 and K 3 = 2.2361). The five localized modes are equally 
spaced by AE^omb = \/k^ - = 1 . 
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where W['0i,'02, V's] is the Wronskian of the s functions 0i, and 

kV[ 0 i,' 02 , •••, V's] nieans that the function should be deleted in the sequence 

(■01, 02 , ■■■,'ips)- The Darboux-Crum theorem states that 


= e^{x) , Hi(j)j{x) = ej(j)j{x) {j = 1,2, ...,N) (21) 


where the new Hamiltonian Hi = —dl + Ui{x) is defined by the potential 

Ui{x) = Uo{x) - 2^1og|>V[</9i,(y92, ...,(pn]{x)\. (22) 

To avoid singularities in the potential Ui{x), the Wronskian W[</9i, 922 , 7 ^ri\{.x) should 

not vanish on the real x axis. In practice the Darboux-Crum theorem ensures that the 
two Hamiltonians Hq and Hi, with potentials Uq{x) and Ui{x) related by Eq.(22), are 
essentially isospectral, apart from the energies ei, 62 , ..., cn that can be added and/or 
deleted from the spectrum of the original Hamiltonian. For example, if Cj {j = 1, 2,..., N) 
do not belong to the point spectrum of Hq and the seed functions (Pj(a:) are unbounded 
as X —>■ ±00, then the associated functions (pj{x) are normalizable eigenstates of Hi and 
thus the effect of the multiple Darboux transformations is to add the energies ej to the 
point spectrum of the original Hamiltonian. A special case is the one corresponding to 
a constant potential for the original Hamiltonian, for which application of the Darboux- 
Crum theorem yields a class of reflectionless potentials of the Kay-Moses kind izu. 
Specifically, let us assume Uq^x) = nf and ei = 0, €2 = — ^ 3 , •••, ev = 

where Ki > H 2 > ■■■ > > 0, and let us take (Pj{x) = exp^Hjx) + (—l)'^Cj exp(— k^x), 

with Cj > 0. The potential Ui{x) of the associated Hamiltonian obtained from the 
Darboux-Crum theorem then belongs to the class of reflectionless Kay-Moses potentials 
EOI. Such potentials can be written rather generally as 


Ui{x) ^ - 2 Alog (det.A) 

where the N x N matrix A is defined by 

-exp|!(K„ + 


(23) 


•^n,m ^n,m T \l 


X 


Kjyi T 

In Eq.(24) are arbitrary real and positive numbers. In particular, with the choice 


(24) 


An 


Kn T ^7, 


m^n 


Kr. 


Kjt 


(26) 


the potential Ui{x) turns out to be symmetric, Ui{—x) = Ui{x). The potential (23) 
sustains N bound states 0j(x) at energies ej = n\ — tA- (j = 1,2, which can be 

found from Eq.(20) or as the solutions of the following linear system of equations |2T] 

0j(x) -|- Aj exp^Kjx) ^ ^ ^ (()u{x) + Aj exp{iKjx) = 0. (26) 

17=1 

The asymptotic behavior of h{x) = W{x) as x —>■ ±cx) turns out to be given by 


h(x) ~ — Ki as X —)■ —00 , h(x) ~ ki as x —>■ cx). 


(27) 
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normalized frequency detuning E 

Figure 3. (Color online) Same as Fig.l, but for a Darboux-Crum transformation 
with N = 5 (ki = i, K 2 = 3.9686, K 3 = 3.8730, and K 4 = 3.7081). The seven localized 
modes are equally spaced by AEcomb = \/uf — K 2 — 0-5- 

This means that the Bragg grating is asymptotically a periodic and inhnitely-long 
grating of strength ki, and the frequencies of the bound modes fall within the stop 
band (—ki, ki) of the grating. 

The Darboux-Crum transformation can provide a useful analytical tool for the 
synthesis of Bragg gratings with target spectral features. It should noted that, as 
compared to other analytical methods of grating synthesis based e.g. on the Gelf’and- 
Levitan-Marchenko inverse scattering method [221 Ea] or on the similar method of 
resonance mode expansion [29], in the Darboux-Crum method the scattering potential 
q{x) does not vanish at a: ± 00 [see Eq.(27)], i.e. the grating is synthesized from an 
underlying uniform and inhnitely-extended grating structure by the introduction of 
defective modes. This is especially useful in the synthesis of DFB cavities that sustain 
a target number of modes. As an example, we apply the Darboux-Crum transformation 
to the synthesis of DFB frequency comb transmission hlters. Integrated optical devices 
such as DFB frequency comb grating hlters can be useful for several applications, for 
example for the minimization of lasing wavelength errors in DFB lasers, for frequency 
selection in multiwavelength excitation and coupling of lasing and nonlinear cavities, and 
for the realization of integrated and spectrally tailored sources of phase-locked lasing 
modes (see, for instance, [HHlEBE21E3] and references therein). To realize an integrated 
DFB cavity sustaining a comb made of {2N — 1) axial modes with frequency spacing 
^^comb = AF'comfe/(2vrT), let US consider a periodic index grating (i.e. 0(a;) = 7r/2) 
and let us apply the previously described synthesis method by selecting the sequence of 
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numbers Ki, ^ 2 , ■■■, Kn according to 

Kj = - [27l{j - l)TAUcomb? (28) 

[j = 1,2,N), where the time scale T is defined by Eq.(4) and where Ki is taken 
arbitrarily, but larger than 27r(iV — 1)TAucomb- The potential Ui{x) and corresponding 
bound states (j)j{x) can be obtained from Eqs. (23-26). The grating profile h{x) = W{x) 
is then retrieved from the expression of the lowest bound state (/>i(x) using Eq.(15). 
A simple expression for h{x) can be given for the case N = 2, i.e. for a DEB cavity 
sustaining three modes: 

— k \ sinh(fi:ix) -|- sinh(Kia;)sech^(K23^) + tanh(K2a^) cosh(Kia;) 
K2 tanh(fi:2a:) sinh(Kia;) — ki cosh(Kia;) 
the frequency spacing being given by 


^ ^comb 


A E(,QjYib 

2'kT 


Ki - nl 


2nT 


(30) 


Figures 1, 2 and 3 show typical examples of grating profiles h{x) for DEB cavities that 
sustain frequency combs made of 3 (Fig.l), 5 (Fig.2) and 7 (Fig.3) modes, as obtained 
from the Darboux-Crum transformations for N = 2, N = 3 and = 4, respectively. 

In a DEB structure of finite length L, the ideal amplitude profile h{x) must 
be truncated after some length far from the ’defective’ region, i.e. h{x) = 0 for 
|a;| < L/{2Z). The effect of grating truncation is that bound states of the potential 
Ui{x) become resonance modes with a finite lifetime (see, for instance, CZD- In a passive 
DEB structure, i.e. in the absence of gain, the resonance modes can be observed as 
narrow transmission peaks embedded in the stop band of the grating, at the normalized 
detunings E = En = {u — ub)T from the Bragg frequency ojb according to Eq.(18). This 
is shown, as an example, in Fig.4(a) for the case of the DEB cavity of Fig.l sustaining 
a comb of 3 modes. The figure shows the transmission spectrum T{E) = \t{E)\‘^ of a 
periodic DEB structure with the index profile h{x) given in Fig. 1(a), taking h{x) = 0 
for |a:| > L/{2Z) = 5. The spectral transmission t{E) versus the normalized frequency 
detuning E has been numerically computed using a standard transfer matrix method 
(see, for instance, [IS]). In Fig.4(a), the three narrow resonance peaks embedded in the 
stop band of the grating correspond to the three resonance states of the comb. Owing 
to grating truncation, other resonance modes appear, however they are high-threshold 
modes and can be disregarded. The resonance modes of the truncated grating can be 
numerically computed by looking for the poles of the spectral transmission t{E) in the 
complex E plane, with E = E — iE^ and Ei > 0 (see, for instance, US]). The real 
part E of the E defines the normalized oscillation frequency detuning of the resonance 
mode, whereas the imaginary part Ei of —E is the corresponding normalized threshold 
gain. For a DEB structure with a uniform gain g per unit length, in physical units the 
threshold gain of each mode is given by 
Ei irAnEi 
^ ^ An 


9th ry 


(31) 
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Figure 4. (Color online) (a) Numerically-computed transmission spectrum {\t{E)\‘^) 
of the DFB filter with index proHle shown in Fig. 1(a) and with h{x) = 0 for |a:| > 5. 
(b) Numerically-computed poles of the spectral transmission t(E) in the complex 
E = E — lEi plane. The real {E) and imaginary {Ei) parts the poles define the 
normalized lasing frequency and lasing threshold of the various DFB cavity modes. 
The inset in (b) is an enlargement of the three low-threshold poles that correspond to 
the three bound states of Fig. 1(c). 


The numerically-computed poles of t{E), i.e. normalized frequency detunings E and 
normalized thresholds E^ of the various DFB cavity modes, are shown in Fig.4(b). Note 
that the three comb modes show a very low lasing threshold [Ei <1 x 10“^; see the inset 
of Fig.4(a)], whereas the other DFB cavity modes have a much higher lasing threshold 
{Ei > 0.015, i.e. more than 20 times higher). To get an idea of the design parameters of 
the DFB structure in real physical units, let us consider a grating realized in a silicon- 
on-insulator (SOI) waveguide [Ml ESI ESI EZ] with Bragg wavelength Aq = 1550 nm. 
Taking uq ~ ~ 3.4, a hrst-order grating period A ~ 228 nm is required [M]- For 

a frequency comb spacing Aucomb = 100 GHz, the refractive index change An and the 
temporal and spatial scales T,Z, as obtained from Eqs.(4) and (30) with ki = 1.2 and 
K 2 = 1, are given by An ~ 5.13 x 10“^, T ~ 1.056 ps and Z ~ 93.2 pm. The full 
grating length is thus L = lOZ ~ 932 pm, the effective coupling constant of the grating 
is ^grating = 7iAnK2/Xo ~ 124.65 cm“^, and the lasing threshold gth of comb modes is 
smaller than 0.0751 cm“^. Such values are realistic for SOI-based waveguide gratings 
fabricated using deep ultraviolet lithography. In particular, a graded-index periodic 
grating with the target prohle as dehned in Fig. 1(a) could be realized by superposition- 
apodized or phase-apodized methods [36], as well as by duty-cycle engineering EH. 

4. Conclusions 

Bragg scattering in grating structures possesses SUSY properties similar to those of 
the one-dimensional Dirac equation. Such an interesting property can be useful for the 
synthesis of DFB hlters and cavities with desired scattering and localization properties 
by application of multiple SUSY (Darboux-Crum) transformations. As an example. 
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we designed DFB structures sustaining with a set of equally-spaced resonances, i.e. 
a frequency comb transmission filter. Further extensions of SUSY methods to Bragg 
gratings could be envisaged. For example, the SUSY structure of Dirac equation in the 
most general case [27] could be considered in the design of DFB cavities with chirped 
gratings. Also SUSY can be applied to underlying non-Hermitian equations, which can 
be useful for the design of DFB cavities that include gain/loss gratings [28]. Further 
extensions include application of second-order (Abraham-Moses) SUSY transformations 

[20] l38] , which might be useful for the design of DFB cavities that sustain bound states 
embedded in the continuum [39] . 
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